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1. Introduction 

The fractional quantum Hall (FQH) effect [1 occurs in two-dimensional interacting 
systems in a strong perpendicular Abelian gauge field at sufficiently low temperatures. 
It manifests as incompressible strongly correlated ground states, at several simple 
fractional values of the filling factor (ratio of particles to flux quanta), with anyonic 
excitations. The FQH effect is usually realized with electrons in a magnetic field, but 
theoretically exists for both bosons and fermions, and both long- and short-ranged 
interactions [2j [3]. Several methods have been proposed for implementing it with 
ultracold atoms [4] [5j |6j [TJ [HJ [9j [TDJ [11] , and some few-atom experiments have been 
done 112]. 

One motivation for studying the FQH effect is that it can support anyons [13] [14] . 
These are particles whose exchange statistics is not the ±1 of bosons or fermions, but 
some other phase e l6> for Abelian anyons, or a unitary matrix for non- Abelian anyons. 
This statistical effect depends only on the topology of the anyons' paths and not on 
the details of their motion. In the case of non- Abelian anyons, the states on which the 
statistical matrix acts are a global property of the pair, indistinguishable if each anyon 
is measured individually, but can be measured by bringing them together. Hence, a 
quantum computer using well separated non- Abelian anyons as qubits and statistical 
matrices as gates would be immune to error from either local decoherence or imperfect 
control of the anyon motion [15] , as long as there were no uncontrolled stray anyons in 
the system. For some types of anyons, including those expected in some FQH states 
P~6j[l3], such a quantum computer would be universal. 

However, errors will occur if uncontrolled stray anyons, such as thermally created 
anyons [13 [18], move around the computational anyons. If the temperature is too high 
compared to the anyon pair creation gap, these stray anyons will make the system 
unusable for computation, and eventually destroy the FQH effect completely. In 
standard FQH systems, the states supporting non- Abelian anyons have small enough 
energy gaps that they are difficult to observe [19] |20] and would be even more difficult 
to use for computation. This provides a motivation to look for modified FQH systems 
with more robust non- Abelian anyons. 

The modification we consider here is to add a uniform SU(2) non-Abelian 
component to the magnetic field, coupled to a two-dimensional internal state space of 
the particles. Such a system can be realized using ultracold atoms in an optical lattice 
[H [HJ [22] [23] [24] [25] [26] , or using electrons in a material with spin-orbit splitting 
[27] [28]. We study this system using lowest Landau level exact diagonalization, for 
bosonic particles and a two-parameter range of SU(2) field strengths. When both 
parameters are of the same order as the Abelian field strength, the system behaves 
similarly to a one internal state FQH system, i.e. incompressible states at filling 
factors v = |, 1, |, . . . with large overlaps with the Read-Rezayi states [29 ], despite 
actually having two internal states. For most parameter values, these have smaller 
energy gaps than in a true one internal state system, but in the cases v = 1, |, it 
appears that some choices of non-Abelian field can slightly increase the gap. When 
both parameters are large, we find no incompressible states at all, while when only 
one is large, we find similar behaviour to when both are zero, i.e. a two internal state 
FQH system. 

In section [2] we introduce our model, starting from a lattice Hamiltonian, then 



taking the continuum limit. In subsection 2.3 we derive the single particle Landau 



levels in the presence of the SU(2) gauge field. Section [3] gives many-particle exact 
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diagonalization results for the density profile, the energy gap, and the overlap with 
some trial states. In particular, we consider the v = 1/2 Laughlin state and the v = 1 
Moore-Read state. We conclude in section HJ 



2. The Model 



We consider bosons in two dimensions, in either continuous space or a lattice that 
can be well approximated by continuous space. We apply a uniform classical 
U(2)=U(1) xSU(2) gauge field, minimally coupled to two internal states of the 
particles. The U(l) part is the Abelian field of standard FQH, while the SU(2) 
non-Abelian part is our new ingredient. The FQH effect also requires an interaction 
between the particles, which we take to be a contact interaction. A possible physical 
implementation of this system is described in [2TJ 0] , using ultracold bosonic atoms in 
an optical lattice with laser induced hopping. 



2.1. Lattice Hamiltonian 

We consider bosons with two internal states on a square lattice, described by the 
Bose-Hubbard Hamiltonian with a classical gauge field [4 , 



u x,y 



+4 >3/ exp ( i / A(x,y')dy' ] a^_i + H.c. 



. y- 

+ Ua i,y a i,y a x,yCLx,y (1) 

Here A(x, y) is the vector potential of the gauge field, cl^ creates a particle in lattice 
site (x,y) (we use the lattice spacing as the unit of length), J is the hopping rate 
and U is the interaction strength. For simplicity, we assume isotropic hopping, and 
state-independent hopping and interaction. 

In the case of a non-Abelian gauge field, the particles have multiple internal (e.g. 
hyperfine) states, and A is a matrix in internal state space as well as a vector field in 
real space, i.e. hopping can change the particle's internal state. We consider the case 
of a uniform U(2)=U(l)x SU(2) gauge field [21], 

A(x, y) = (ao-y, j3cr x + 2tt§x) 

-ia \ / 2tt$x (3 



ia )'\ $ 2tt^x 1 1 ' ( 2) 
where the are the Pauli matrices. The Abelian part of this field is parameterized by 
<£, and the non-Abelian part by a,/3. As we see below, the Abelian part is necessary 
to have an FQH effect; a pure SU(2) field [30] does not produce Landau levels. 

For this choice of field, A x is independent of x and A y of y, so the integrals in ([T]) 
simplify to e lAx , e lAy , giving 

H = ^ — J [a^y (cos a + ia y sma)a x -i yy 

x,y 

+e 27ri ^a^(cos/3 + ia x smf3)a x , y -i + H.c] 

^^ a x J y a x J y a x,y a x y y- (3) 
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2.2. Continuum approximation 



When the lattice spacing is small compared to the magnetic length (defined below) 
and the typical interparticle spacing, we can approximate the site operators a XiV by 
a continuous field ip(x,y). Taylor expanding ip{x — l,y) and similar terms in ([3| 
about (x,y), and discarding the (dynamically irrelevant) constant —AJip^ip, gives the 
approximate Hamiltonian 

2 



jdxdyJ^ ( 



.d_ 

dx 



-i— + f3a x + 2tt<$>x 
dy 



ip -\ — Uip^ip^ipip, 



(4) 



where all occurrences of ip are taken at position (x,y). This Hamiltonian has a 
magnetic length / = l/v / 27r<i>, which sets the length scale of the quantum Hall physics. 
Hence, the above continuum approximation is valid in the limit I ^> 1. For convenience 
we also define the dimensionless parameters a = a/y/$> and b = /3/y/$>. Within this 
continuum limit, varying at fixed a, b changes only the overall length and energy 
scales. We can assume without loss of generality that a, 6, $ > 0, as a spin space 
rotation by a x (a y ) changes the sign of a (6), while a real space reflection in the y axis 
changes the signs of a and <£>. 

Expanding out the brackets and discarding the dynamically irrelevant constants 
a x — a y ~ ^ we fi n d that the non-interacting part of Q is equal to the Hamiltonian of 
2D electrons with spin-orbit coupling [28], with the coordinate axes rotated by 45° (in 
both real and spin space). The SU(2) terms correspond to the spin-orbit interaction, 
which is pure Rashba if a = — /?, pure Dresselhaus if a = /?, and a combination of the 
two for all other ratios. An optical lattice loaded with fermionic atoms can hence be 
used as a quantum simulator of spin-orbit coupled electrons [4 . However, we consider 
only the bosonic case. 



2.3. Single particle degeneracy: Landau levels 

In this subsection we consider the behaviour of a single particle subject to the 
Hamiltonian Q. Since H c is independent of its single particle eigenstates are 
plane waves in the y direction, ip ex e lky . Substituting this into we find that the 
eigenstates have the form ip = e lky ty n (x + k/27r<&), where \I/ n and the energy E n satisfy 
the equation 

E n V n (X) = J[(-id/dX + ao- y f 

+(/3a, + 2^X) 2 ]^ n (X). (5) 

Like the Abelian case, the energy is hence independent of so each single particle 
energy level n (called a Landau level) has a macroscopic degeneracy of $ states per 
unit area (one per Abelian flux quantum). In the many-particle case, this degeneracy 
is broken by the contact interaction, so a weak interaction can produce strongly 
correlated states. As in standard FQH, we define the filling factor v as the number of 
particles divided by the number of states in a Landau level, v = p/§ where p is the 
2D density. 

In the Abelian case (a = f) = 0), E n = 2, 2, 6, 6, 10, 10, . . . x ttJQ, the additional 
2-fold degeneracy coming from interchange of internal states. In contrast, in the non- 
Abelian case, the Landau level energies E n are in general all distinct and follow no 
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Figure 1. Lowest few Landau level energies E n , and symmetries H — (blue/dark) 
or — h (green/ light), for a = b (thin lines) and a = b ± 1 (thick lines). 

simple pattern. The E n can be found analytically [28], but this result is complicated, 
and gives the corresponding wavefunctions \I/ n only as infinite series; we instead solve 
([5| numerically. Full analytic solutions exist in the special cases a = or f3 = (2-fold 
degenerate [27]), or a = ±(3 (non-degenerate [3T| 132]). 

Figure [T] plots the lowest few Landau level energies E n along two diagonal lines 
in (a, b) space, while figure [3^t plots the difference between the lowest two energies 
against a, b. Figure [2] plots the lowest Landau level wavefunction \I/o at four (a, b) 
settings. Because ([5| has the symmetry ^ n (x) —> o- z ^! n (—x), in each Landau level 
either the first internal state component is symmetric about x = — k/2ir<& and the 
second antisymmetric (H — symmetry, figure [2^i,b,c), or vice versa ( — h, figure |2]i). 
Levels of the same symmetry avoided cross, touching only in conical intersections on 
the a = ±6 diagonals, but levels of opposite symmetry cross freely, on lines in (a, b) 
space. For small positive a, 6, the lowest Landau level is H — and the second — h, but 
at larger a, b these cross infinitely many times [32 , causing a discontinuous change in 
the wavefunction each time. 

The exactly solvable cases a = and b = can be used as perturbative 
approximations for the cases when a or b is small but not zero [27]. When a = 0, 
the two degenerate ground states at a given k are Gaussians with shifted centres, 
ifj(x,y) = exp{i% — [x + (k ± /3)/27r<£] 2 /2/ 2 }, in the a x = ±1 internal states (equal 
superpositions with relative phase ±1). When a <C 6, the ground state is close to the 
symmetric combination of these (figure [SJd) , and hence has two peaks approximately 
fS/irQ = bl/y/n apart (figure [3]i) . When 6 = 0, the ground states have a Gaussian 
amplitude but an extra phase, ip(x,y) = exp[i(% ± ax) — (x + k/2it<&) 2 /2l 2 ] with 
internal states <j y = ±1. The symmetric combination, and approximate a ^> b ground 
state (figure^), has a Gaussian total amplitude but a position dependent internal 
state. 

However, this apparent asymmetry between a and b arises entirely from our choice 
to have states extended along y and localized in x, and to use a gauge that simplifies 
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Figure 2. Single particle lowest Landau level wavefunctions ^o(x), for non- 
Abelian field strengths (a, b) = (a) (1,1), (b) (1,5), (c) (5,1), (d) (5,5). Blue 
and green are the two internal states ^o(^t)> ^o(x,i), black their rms total 
amplitude y / |^o(^,t)| 2 + |^o(^,^)| 2 , thin grey the Abelian-field ground state for 
comparison. 



such states; we could instead have made the Abelian part of A proportional to y (by 
a gauge transformation) and had states ip ex e lkx , which would interchange the roles 
of a and b. 

3. Numerical results 

We now consider the multi-particle case. We assume weak interaction, by which we 
mean that the interaction energy per particle (~ U$> at v ~ 1) is small compared to 
the Landau level spacing (~ i.e. [/ « J. We can hence assume that all the 

particles are in the lowest Landau level, and use a Fock basis over the single particle 
eigenstates \k). Since all momenta k within a Landau level are degenerate, the only 
non-trivial term is the interaction term, with the matrix elements 

/juout jLOutl TT I i.in iin\ _ T]Z(\M i L,out Z~out\ 

i ^2 | -°int 1^1 5 ^2 / ~ u °\ hlJ l * ^2 ^1 ^2 / 

x J dx%(x + k^^2n$)%(x + k% ut /27r$) 

x^ (x + k[ n /27r<S>)^ (x + K/2tt^). (6) 

We adopt periodic boundary conditions (a torus geometry) [33j[34], which makes the 
number of allowed k finite, and use the full 2D Haldane symmetry [34]. The interaction 
Hamiltonian then becomes a finite matrix, which we diagonalize to find the ground 
and lowest few excited states. This method is exact within the lowest Landau level, 
and in particular can describe the strong correlations typical of FQH states, but is 
limited to small particle numbers by the exponentially growing basis size. 
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Figure 3. Single particle state properties: (a) Energy difference (E\ — Eq)/Ai:J& 
between the lowest two Landau levels, and symmetry (H — or — h) of the lowest 
level, (b) Overlap {^ . Abelian |^o) of the lowest Landau level wavefunction with 
the Abelian lowest Landau level, (c) RMS spread (standard deviation) of the 
lowest Landau level wavefunction, ((^o|^ 2 |^o)) 1//2 A- The small dark features 
on the right hand edge are numerical artefacts, caused by the gap between the 
two lowest Landau levels being too small to resolve, (d) Distance between the 
two highest maxima of (equal height by symmetry, see figure |2j),d), or zero if 
there is a single central maximum. 
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Figure 4. Filling factor (zeros offset for clarity) against chemical potential, found 
by exact diagonalization on a square torus, (a) Abelian field (single internal state), 
(b-d) non- Abelian fields (a, b) = (b) (2,2), (c) (5,5), (d) (5,1). System sizes 6 (red), 
8 (yellow), 10 (green), 12 (cyan), 16 (magenta) flux quanta; due to computational 
limitations, the larger sizes can only be used at small filling factors. Horizontal 
lines mark the Read-Rezayi densities u = 1/2, 1, 3/2, .... 
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Figure 5. Filling factor (lines mark plateau edges, Read-Rezayi series v = p/2 
labelled in black text, composite fermion series v = p/(p db 1) in blue, others 
in magenta), and projected overlap with the corresponding Abelian-field state 
(shading), against chemical potential and a. Overlap is with the lowest Abelian 
state of the same Haldane momentum K, projected into the non-Abelian lowest 
Landau level; green crosses indicate that this is not the same K as the Abelian 
ground state, (a) a = 6, (b) 6=1. System sizes 12 {y < 1, cyan lines), 8 
(1 < v < 2, yellow), 6 {y > 2, red) flux quanta, on a square torus. Dark red 
triangles indicate densities skipped because the exact diagonalization failed to 
converge. 



3.1. Density profiles 

We perform exact numerical diagonalization for different numbers of particles (at 
a fixed number of Abelian flux quanta, i.e. number of single particle states per 
Landau level) to obtain the ground state energy as a function of density. We then 
convert this to density as a function of chemical potential [35] . This approach makes 
incompressible states visible as plateaus in the density. Furthermore, in the local 
density approximation (valid for weak traps) the local chemical potential equals the 
global chemical potential minus the trap potential [35] , so these results can be directly 
compared with experiment if the in-trap density profile can be measured. 

At moderate non-Abelian field strengths (a, b < 4), the density profile (figure [4]d, 
lower halves of figures|5|i,b) resembles the single internal state Abelian-field case (figure 
[4^i and [35]). Furthermore, most of the ground states have large (> 90%) overlaps with 
the projections of their Abelian-field equivalents onto the new lowest Landau level. 
This makes sense, as though we have two internal states, the lowest Landau level is 
a single level, not a degenerate pair. The main plateaus are the Read-Rezayi series 
v = p/2 (p integer), clearly present for v = 1/2, 1, . . . , 3 and possibly present all the 
way up to our calculation's limit of v = 4. (Because the number of particles must be 
an integer, finite system density profiles are stepped even in compressible phases, so 
narrow incompressible plateaus cannot be reliably distinguished from regions of low 
but non-zero compressibility.) Of the composite fermion series v = p/(p± 1), we see 
hints of v = 2/3,3/4,4/5,7/6,6/5,5/4 (but not 5/6 or 4/3). We also see v = 7/4, 
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9/4 and 11/4, but do not have any theory as to what these states are (the previous 
Abelian- field study [35] used 6 flux quanta, at which these v are not possible). They 
may be finite size artefacts, as the 5/4 plateau narrows greatly on increasing the 
system size from 8 to 12 states. Due to computational limitations, we were unable to 
do this test for the higher densities. There are two main differences from the Abelian 
case at these values of a, b. First, all plateaus are shifted towards higher chemical 
potential. This would cause a trapped system (with fixed particle number) to slightly 
expand. Second, the Laughlin state (y = 1/2) no longer has exactly zero interaction 
energy, and possible new plateaus appear below it at v — 1/4, 1/3, though again these 
are too narrow to be sure they are real. 

If the non-Abelian field strength is further increased along the a = b diagonal 
(figure [i]:, upper half of figure^), all the plateaus disappear suddenly, in a first order 
phase transition. As this happens at the same point as the single particle Landau 
levels cross (figure [I]), it is probably caused by this sudden change of single particle 
states. 

If a is increased while b is kept constant (figure [4]i, upper half of figure [5]o) or vice 
versa, most of the plateaus disappear gradually, though at some system sizes a few 
plateaus survive or even grow (y — 1/3, 2/3, 1 for 12 flux quanta), then at even larger 
a the normal plateaus reappear. This process becomes slower and occurs at a larger a 
as the system size is increased, with the minimum plateau width occuring very close to 
the point where the single particle peak separation equals half the torus circumference 
(a = 4.1,4.6,5.5 for 6,8,12 particles and 6=1, figure |3jl). We hence believe the 
revival, and possibly some of the initial decrease, to be a finite size artefact, caused 
by single particle states wrapping right round the torus. In the infinite size limit, it is 
not clear whether this will become a second order phase transition to a gapless phase 
at finite a (as suggested by simple 1/N extrapolation, figures [6pb, [8pd), or a decay 
that tends to zero only in the a —> 00 limit. 

However, as the gap between the lowest two Landau levels is exponentially 
decreasing in this region (~ 5% of its peak value at the point of minimum plateau 
width, figure [3|i), the weak interaction assumption must eventually fail, allowing the 
second Landau level to be populated. In the a > 6 or 6 > a limit, the two Landau 
levels become degenerate, and by taking superpositions of them we can recover two 
sets of single particle states with Gaussian amplitudes. The two states centred at 
x = xq are ijj(x,y) = exp(— (x — xo) 2 /2l 2 — 2ixoy/l 2 =p iax) (a ^> 6, internal states 
(j y = ±1) or ijj(x,y) = exp( — [x — xq) 2 /21 2 — 2\x§yjl 2 =p i/3y) (a <C 6, internal states 
a x = ±1). Because the interaction does not change the internal state, the extra 
phase factors cancel out in the interaction Hamiltonian ^ ijj^ i/ji/j , so in this limit the 
Hamiltonian matrix is identical to that at a = b = 0, i.e. Abelian-field FQH with two 
internal states. This is known to have a v = 2/3 221 state [2], and is conjectured to 
have a v = 2p/3 NASS series [36] . 

The weak interaction assumption would also fail close to Landau level crossings 
[32] , but here the two degenerate single particle states are not close to Gaussian. We 
hence do not know whether these parameter values would be in a two internal state 
FQH phase, or a gapless phase similar to the one internal state gapless phase above 
the first such crossing. 
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Figure 6. v = 1/2, (A) N = 6 particles, (B) N = 8 particles, (C) 1/N 
extrapolation to infinitely many particles, based on data for 4-8 particles, (a) 
Overlap of the numerically calculated lowest K = (0, 0) state with the projection 
of the exact Laughlin state onto the lowest Landau level, (b) Energy gap between 
this state and the first excited state (of any K). (c) Energy gap between the 
excitation branch and the continuum (as shown in figure [T^l) , smallest gap of 
K = (N/2, 0), (0, N/2), (N/2, N/2) (K scaled to take integer values 0, . . . , N- 1). 
The small features near (a, b) = (4,8), (8,4) are numerical artefacts. 



3.2. v — 1/2: Laughlin state 

We now consider the ^ = 1/2,1 states in more detail. The bosonic Laughlin state, 

Hz 1 ,...,z N ) = H(z i -z j fe-^ 2 ^ (7) 

where zi = Xi + ijji is the position of particle i as a complex number, has exactly zero 
interaction energy as it is zero whenever two particles come together. In Abelian field, 
it is also a lowest Landau level state, so is the exact ground state at v = 1/2. It is 
a non-degenerate state (except for the 2-fold centre of mass degeneracy of all states 
at half-integer filling [34]) with zero Haldane momentum, with a gap to a well-defined 
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Figure 7. Lowest few energies found by exact diagonalization, plotted against 
\K\. (a) v = 1/2 (8 particles, 16 flux quanta), (a, b) = (2,2). (b,c) v = 1 
(12 particles, 12 flux quanta), (a, 6) = (b)(2,2), (c)(2,l). (d) i/ = 3/2 (15 
particles, 10 flux quanta), (a, 6) = (2,2). Colour denotes projected overlap with 
the corresponding Abelian-field state, from none (red) to full (black). The green 
circles mark the ground state (s) predicted by Read-Rezayi theory. 



branch of excitations [2|. 

In non-Abelian field, the Laughlin state is no longer a lowest Landau level state, 
but it can be projected onto the new lowest Landau level and used as a trial state, 
which has large overlap with the ground state for nearly all a, b (figure^). Its 
interaction energy is no longer zero after this projection, but for a, b < 4 the energy 
spectrum remains qualitatively similar (figure [T^i) , retaining the gap (figure ^p) and 
well defined excitation branch (figure [6^) . 

At stronger fields, if a w b the gap suddenly disappears, while for a ^> b or 
a <C b it continuously decreases to near zero then increases again. As discussed in 
subsection |3.1| for general z/, the first is probably a first order phase transition driven 
by the sudden change of single particle state, while in the second case, it is not clear 
whether or not the gap actually reaches zero at finite field strength, while the revival 
is probably a finite size artefact caused by single particle states wrapping round the 
torus. 

Our method does not reveal the excitation statistics. However, the excitations of 
any non-degenerate gapped state at non-integer filling have fractional particle number 
by Laughlin's gauge argument [l], which is still valid in the presence of the non- 
Abelian field, and hence fractional statistics [37] . This does not distinguish Abelian 
from non-Abelian statistics, but given the lack of a visible phase transition we expect 
the same statistics as the Laughlin state, i.e. Abelian = tt/2 (semions) [38] . 

3.3. v = 1: Moore-Read state 

The pth Read-Rezayi state for bosons is defined as the v = p/2 state which vanishes 
when any p + 1 particles come together; it is unique on a disc but p+ 1-fold degenerate 
(including the twofold centre of mass degeneracy if p is odd) on a torus [40] . It is equal 
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Figure 8. v = 1, (A) N = 8 particles, (B) N = 12 particles, (C) 1/N 
extrapolation to infinitely many particles, based on data for 6, 8, 10 and 12 
particles. Due to strong size dependency, this extrapolation is not reliable, 
(a) Haldane momenta K below the largest energy gap in the lowest 10 states. 
White: Moore-Read (K = (N/2, 0), (0, N/2), (N/2, N/2)). Cyan (pale): Smectic- 
Y (two or more K with equal K y and evenly spaced K x , suggesting stripe order 
parallel to y [39]). Red (medium): Smectic-X (two or more K with equal K x 
and evenly spaced K y ). Blue (dark): Single state (usually, but not always, 
K = (N/2,0) or (0,N/2)). Black: any other combination. (b,c) Overlap of 
the numerically calculated lowest K = (N/2, N/2) (b) or K = (0,N/2) (c) 
state with the projection of the exact Moore-Read state onto the lowest Landau 
level. For the third component (K = (N/2,0)), reflect (c) in the a = b 
diagonal, (d) Energy gap between the highest of the three Moore-Read states 
(K = (N/2, 0), (0, N/2), (N/2, N/2)) and the next lowest state (as shown in figure 
[7p); negative (truncated to zero in the figure) if these three states are not the 
lowest three, (e) Energy splitting between the highest and lowest of these 3. 
Lines in (d,e) are energy level crossings of these 3 states, cyan if the ground state 
is involved, red if not. 



to the symmetrized product of p Laughlin states [41] . It is not an exact eigenstate for 
p 7^ 1, but in Abelian field has been found to be a good trial state for small p [29] . 

The p = 2 (v — 1) case is also called the Moore- Read or Pfaffian state. 
Its excitations are Ising anyons, whose non-Abelian statistics have been tentatively 
detected in a related fermionic state [20]. These are not computationally universal on 
their own [42], but can be made so with relatively noisy non-topological operations 
[42] [43] . On a torus the Moore- Read state is 3 states at Haldane momenta 
K = (N/2, 0),(0,JV/2), (N/2, N/2) [29] (where K is scaled to take integer values 
0, . . . , N — 1); for a system to be in the Moore-Read phase, in the infinite size limit 
these 3 states must be topologically degenerate [44] and have an energy gap between 
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them and all other states. However, the Moore-Read state has a smaller energy gap 
and a longer correlation length than the Laughlin state, and hence is more sensitive 
to finite size error, which can split the degeneracy. At the sizes we can access (up to 
12 particles) it is not well converged, as shown by the strong system size dependence 
in figure [8] 

We find a positive energy gap (i.e. the three lowest energy states are at the right 
three K for a Moore- Read state, figure [8]i) and large overlaps between the exact 
ground state and the projected (into the lowest Landau level) Moore-Read state at all 
three K (figure [8)3, c) over roughly the same region as the Laughlin state, i.e. a, b < 4 
plus probably- artefact revivals when a or b is large. However, the gap and degeneracy 
splitting (figure [8^) vary significantly within this region, with 2-3 "good" stripes (large 
gap, small splitting, energy spectrum shown in figure [TJd), separated by "bad" stripes 
(small gap, large splitting, figure [7J3). Simple 1/N extrapolation suggests that only the 
"good" stripes are in the Moore- Read phase (zero splitting in the infinite size limit, 
figure [8pe), but the strong system size dependence makes this result unreliable, and 
it is a plausible alternative that the entire a, b < 4 region is in the Moore- Read phase. 
The stripes arise from crossings of the three Moore-Read state energies (marked by 
lines in figure [8]i,e); near triple crossings, giving a greatly reduced splitting, seem to 
be more common than would be expected by chance. 

At stronger fields, the Moore- Read state is replaced by smectic (stripe ordered), 
non-degenerate gapped, or gapless states (figure [8^1). Which of these occur in what 
regions of a, b parameter space is strongly system size dependent, and we hence do 
not know which, if any, of these survive in the infinite size limit. As discussed above 
for general at these fields the width of a single particle eigenstate is a large fraction 
(> |) of the the torus circumference, and strong finite size artefacts are hence likely. 

We did not perform a similarly detailed study of the v = 3/2 state (which if Read- 
Rezayi supports universal topological quantum computation [16]), as at the sizes we 
can reach, the two K = states that should be its degenerate ground states are 
mostly not even the lowest two states (figure [7]i); this is consistent with previous 
Abelian-field results [45]. However, we note that it too can have its splitting reduced 
and its "gap" made less negative by some choices of non-Abelian field, and that its 
components sometimes avoided cross (one such crossing causes the rapid decrease of 
overlap at a « 2 in figure^), suggesting similar behaviour to v = 1. 

3.4- State- dependent atomic interactions 

The above assumes equal inter- and intraspecies interaction, — Un = U^. For 
hyperfine states of real atoms, these are usually nearly but not exactly equal, and 
they can be far from equal near a Feshbach resonance. Also, an alternative method of 
generating the U(2) field [25] uses the two pseudo-internal states (momentum valleys) 
that occur at ~ 1/2, which have a fixed 1:1:2 interaction ratio [46], instead of 
true internal states. 

We tested the effect of asymmetric interactions by repeating some of the above 
calculations for U n : U u : U n = 1.024 : 0.973 : 1 (the \F,m F ) = |1, -1) , |2, 1) states 
of 87 Rb [2 ), 1 : 1 : (i.e. no interspecies interaction) and 1:1:2. The small 
deviation from 1 had no visible effect; since there is only one lowest Landau level, not 
a degenerate pair, there is nothing to phase separate. The larger increase (decrease) 
of interspecies interaction increased (decreased) the parameter range over which FQH 
was seen and its energy gap, but did not alter its qualitative features. 
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4. Conclusions 

We numerically studied the fractional quantum Hall effect in bosons in a U(2) gauge 
field. For moderate non-Abelian field strengths, we find that it behaves similarly to 
a one internal state Abelian-fleld quantum Hall system, despite actually having two 
internal states. Within this regime, the energy gap of the non-Abelian anyon states 
v = 1,3/2 is strongly dependent on the non-Abelian field parameters, with some 
settings giving larger gaps than in the true one internal state system. However, as the 
gap is also strongly dependent on the system size, it is not clear whether this effect 
persists in the infinite size limit. 

In strongly asymmetric non-Abelian fields, the system reverts to behaving as a two 
internal state quantum Hall system (i.e. like it would without any non-Abelian field), 
while in strong near- symmetric non-Abelian fields, it does not exhibit the fractional 
quantum Hall effect at all. 
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